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This set of exercises relates to Book 3 of SM358 and can be used for further practice and for revision purposes. To 
gain maximum benefit, it is important to attempt each question for yourself before looking up the solution. 


Exercises 
Topic 1 — Spherical harmonics 
Exercise 1.1 (a) Show that the function 
f(@) = Acos@, 
where A is a constant, is simultaneously an eigenfunction of 


| Ope tO 
_ E sin 65 + Sin? 6 O@ 


and 
ES oO 
L, = = oy 
Find the corresponding eigenvalue in each case. 
(b) Find a value of A that ensures that f (6) is a normalized spherical harmonic, Yj,,,(0,@). What are the values 
of l and m in this case? 
You may use the standard integral 


f cos? 0 sin 8 dé = a 
0 3 
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Exercise 1.2 At time t = 0, a system is described by the normalized wave function 


U(r, 0) = f(r) [2%11(8, $) Ez 4¥i0(0, d)| 
(a) Show that the angular part of the wave function is normalized. 


(b) If the square of the orbital angular momentum is measured at t = 0 in the state U(r, 0), what is the 
probability of getting the value 2h?? 


(c) Ifthe z-component of the orbital angular momentum is measured at t = 0 in the state U(r, 0), what is the 
probability of getting the value A? 


Topic 2 — Total angular momentum 


Exercise 2.1 Show explicitly that 


IA) = Nol 1) — 2 Yil) 
IB) = 2 Molt) + yira 


are both normalized and are orthogonal to one another. 


and 


Exercise 2.2 Show that 


IB) = YZY ol t) + VE Mal) 


is an eigenvector of J; and find the corresponding eigenvalue. 


Exercise 2.3 (a) The identity 

Tem a a 
can be rearranged to give 

aasa IR +. z2 

L-S=5( -L -S ). 
Use this result to show that the kets 

ljm) =|1, 5,5) and |j m) = |1, 3, 5), 


which represent states of a spin-4 particle, are eigenvectors of the operator L - S for such a particle, and find the 
corresponding eigenvalues. 


(b) Why, in the physics of atomic hydrogen, is it significant that these eigenvectors have different eigenvalues of 
L.-S? 


Topic 3 — Hydrogen atom 


Exercise 3.1 Allowing for different values of the spin magnetic quantum number ms, determine how many 
degenerate states in the Coulomb model of a hydrogen atom have energy — Eg /16, where Ep is the Rydberg 
energy. Explain how you obtained your answer from the restrictions that apply to the quantum numbers l, m and 
Ms. 


Exercise 3.2 In the Coulomb model of a hydrogen atom, the energy eigenfunctions are written as Ynim(r, 9, @), 
where n is the principal quantum number, / is the orbital angular momentum quantum number and m is the 
magnetic quantum number. State the mathematical and physical reasons why the only allowed values of m, l and 
n are integers. 


Exercise 3.3 In the Coulomb model of a hydrogen atom, radial functions with l = n — 1 take the form 
Ryy(r) = Cr? te77/7 20, 


where C is a normalization constant and ag is the Bohr radius. Show that the radial probability density for a state 
with 1 = n — 1 has its maximum value at r = n?a. Use your answer to estimate the typical size of a hydrogen 
atom with l = 99 and n = 100 and calculate how much energy would be needed to ionize the atom in this state. 
The Bohr radius is ag = 5.29 x 1071! m and the Rydberg energy is Eg = 13.6 eV. 


Exercise 3.4 The energy eigenfunction in the 2s state of a hydrogen atom takes the form 
1 
VAr 


where the radial function is 


1 r 
= _ —r/2a 
pate (2ag)3/2 (2 z) ° ” 


and ag is the Bohr radius. Show that the expectation value of r is 6ao. 


#2,0,0(7, 8, ¢) = Rəo(r), 


You may use the standard integral 
OO 
f re /4 dr = n! att, 
0 
for any positive constant a and non-negative integer n. 


Exercise 3.5 (a) Find the expectation value of 1/r? in the hydrogen atom state described by the normalized 
energy eigenfunction 


1 fey 
#2,1,0(7, 9, 9) = re (=) eo cos 6, 

where ag is the Bohr radius. 

You may use the standard integrals 


T » [0,6] 
f cos? 8 sin ô dé = 3 and f rre"/4 dr =n! ar: 
0 0 


where a is a positive constant and n is a non-negative integer. 


(b) Can you suggest a physical context in which evaluating the expectation value of 1/r? in a hydrogen atom 
energy eigenstate would be useful? 


Exercise 3.6 This question concerns the hydrogen atom energy eigenstate 
WnimlT, 0, p) = Feil?) Vil? ¢) E Rri(r)O1m(0)e"?. 


(a) Show that the expectation and uncertainty of r are independent of the magnetic quantum number m in the 
state nim- 


(b) Show that the probability density at any given point is independent of the azimuthal angle ¢ in the state Ynim- 


(c) Show that the radial probability density at a given radial coordinate is independent of m, @ and ¢ in the state 


Wnim- 


Exercise 3.7 This question concerns the ground state of a hydrogen atom, for which the normalized energy 
eigenfunction is 


1 1/2 
W1,0,0(T, 0, $) = (=) eo 7/20, 
TA 


0 
where ag is the Bohr radius. Given that the operator for the square of the momentum of a particle is 


AD 
a 2_ 2| 190 20 L 
p = ev? =e | r2 Or ” Or aoe i 


show that 
2ao 


R2 
P'v10,0(7, 9, $) = a2 (= = 1) #1,0,0(7, 0, $), 
0 


r 
and hence find the expectation value of p? in the ground state of a hydrogen atom. 


You may use the standard integral 
i 1 
1 re t/* dr = nla"t!, 

0 


for any positive constant, a and non-negative integer n. 


Topic 4 — Perturbation theory 


Exercise 4.1 A particle of mass m is in the one-dimensional infinite potential energy well 


Vía) A|z| for —L/2 < x < L/2 
a e = 
ee) elsewhere, 


where A is a small constant. 


Use first-order perturbation theory to estimate the lowest two energy eigenvalues for this system, taking as the 
unperturbed system the particle in an infinite square well with walls at x = +L/2. 


You may use the integrals 


T/2 2_4 T 2 
i ucos? udu = — and f usin? udu = —. 
i 16 A 4 


Exercise 4.2 A particle of mass m is in the one-dimensional infinite potential energy well 
À for 0 < |z| < w/2 
V(r) = 40 for w/2 < |z| < L/2 
ee) elsewhere, 


where À is a small constant and 0 < w < L. Use first-order perturbation theory to estimate the ground-state 
energy of this system. 


You may use the standard indefinite integral 


1 
fos udu = 5 + Z sin(2u). 


Exercise 4.3 For the purposes of this question, suppose that the electrostatic potential energy of two point 
charges e and —e, separated by a distance r, is given by 


2 —r/X 
e° 1l+e 
V(r) = -——_——_, 
ATED i 
where A is a positive constant that is small compared to the Bohr radius, ag. What consequences would such a 
deviation from Coulomb’s law have for the ground-state energy of a hydrogen atom? 


You can model the electron and proton in a hydrogen atom as point charges, and take the ground-state energy 
eigenfunction in the Coulomb model to be 


141 
¥1,0,0(7, 9, ¢) = (=) es, 
TA 


0 
You may use the standard integral 


oo ! 
n: 
re dr , 
0 a” 1 


where q is a positive constant and n is a non-negative integer. 
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Exercise 4.4 The relativistic perturbation to the kinetic energy of a particle of mass m can be approximated by 
the expression 
2... (py 
ôH = — Sa: 
Find the corresponding first-order correction to the ground-state energy of a hydrogen atom. You may use the 
following facts: 


e The unperturbed ground-state energy eigenfunction is 
1 \12 
Yı o,0(r, 0, ¢) = ( =) eet, 
where ao is the Bohr radius. 


e Because p” is Hermitian, we have (4)|p*|¢)) = (P Yp y). 


2 
x perrea = (280) ef 
0 


You may use the standard integral 
Co 
f re"/4 dr = n!a”t! 
0 
for any positive constant, a and non-negative integer n. 


The following exercise asks you to derive a useful result that is not discussed elsewhere in the course. It is 
included here because some of the methods used in the derivation are similar to those used to obtain the formula 
for the first-order correction in perturbation theory. 


Exercise 4.5 Suppose that a Hamiltonian operator H à depends on a parameter A, which may take different 
values in different systems. For example, the Hamiltonian operator for a harmonic oscillator 


he d? 1 


H =-—— + -Cz? 
tide 2 ay 


contains the constants m and C’, which can be regarded as parameters. Given the Hamiltonian operator fi à» the 
corresponding energy eigenvalues and normalized eigenvectors also depend on the parameter À, and may be 
denoted by En, à and |n, à), so we have 


Hy|Un,a) = Byala): 


In this question, you are asked to show that 


OE. dH 
a= (Yna AÀ a Vn x): 
a result known as the a n theorem. Show first that 
OEn 1 OVn, ofi Yn 
= (> = |, Vn ») + (Yna tna) + (na A, | 5 
Then use the fact that H à is Hermitian, and the normalization condition (w,,y|¢n,,) = 1 to show that 
Ow Ow 
( alfi Úna) + (Yna fix AN =0, 


and hence complete the proof. 


Having derived the Hellmann—Feynman theorem, it may be interesting to see it in use. 


Exercise 4.6 (a) The Hamiltonian operator for a harmonic oscillator is 


where C is the force constant and m is the mass. This has energy eigenvalues 


Ey = (n+p E 


with corresponding eigenfunctions w,,(2), where n is the quantum number. Treating C as a parameter, use the 
Hellmann—Feynman theorem to show that the expectation value of the potential energy in the eigenstate ~,, is 
equal to half the total energy En. 


(b) The Hamiltonian operator in the Coulomb model of a hydrogen atom is 


a O a 
2u Aneor’ 
where e is the electron charge and u is the reduced mass. This has energy eigenvalues 
ef u 1 


[a E 
7 (4reo)? 2h? n?’ 


with corresponding eigenfunctions Vnim(r, 0, ¢). Treating the electron charge e as a parameter, use the 
Hellmann—Feynman theorem to show that the expectation value of the potential energy in the eigenstate Ynim is 
equal to twice the total energy En. 


Topic 5 — Variational method 


Exercise 5.1 A one-dimensional system consists of a particle of mass m, subject to the finite square well 


potential energy function 
— Vo forO<a< JL 
V(x) = 


0 elsewhere. 


—ax? /2 


(a) Use the trial function @ = e , where a is a positive adjustable parameter, to show that 


A 2 vaL 
p, = HRl) Vo [e ay, 


(iloi) Am VT Jo 


is the Hamiltonian operator of the system. 


You may use the standard integrals 


[0,6] OO 
ee dr =,/Z and ae dr = v. 
= Qa = 2a3/2 


where q is a positive constant. 


(b) Show that, when a is very small and positive, Æ, < 0, and hence show that the ground-state energy 
eigenfunction in any finite square well has an energy below the top of the well, no matter how shallow or narrow 
the well may be. 


You may use the fact that 


zx 
f oe du~g when z <1. 
0 


Exercise 5.2 A particle of mass m, in three dimensions, has the Hamiltonian operator 


F h? 2 1 2m2 
H= zao + gmr 5 


where wg is a positive constant and r is the particle’s radial coordinate in spherical coordinates. Obtain a 
variational estimate for the ground-state energy of this system, using the trial function 


g(r) =e", 
where a is a positive adjustable parameter. 


You may use the fact that, when V? acts on any function f(r) that is independent of 0 and ¢, it gives 
1d df 
2 — ee: foe 
Ya)= r? dr (« 2) l 


You may use the standard integral 
oO 
f ree t/* dr = nla” t., 
0 
for any positive constant a and non-negative integer n. 


Exercise 5.3 Repeat the calculation for the Hamiltonian operator given in Exercise 5.2, but using the trial 
function 


plr) = ee, 


You may use the standard integrals 


OO OO 
2 /q2 T 2 a2 3/7 
ret /@ dr = wz a? and ret (© dr = va a”, 
0 4 0 8 
where a is a positive constant. 


Topic 6 — Hydrogen-like systems 


Exercise 6.1 Calculate the scaled Bohr radius (in m) and the scaled Rydberg energy (in eV) for a muonic atom 
of “He, which can be considered to be a bound state of a muon (of charge —e and mass 207me) and a helium-4 
nucleus (of charge 2e and mass 7360me), where e and Me are the magnitude of the charge and the mass of an 
electron. 


You may take the Rydberg energy and Bohr radius to be Ep = 13.6 eV and ay = 5.29 x 1071! m. You may also 
take the reduced mass of a hydrogen atom to be the same as the mass of an electron. 


Exercise 6.2 Find the scaled Bohr radius and the energy of the ground state of positronium, a bound state 

of an electron and its antiparticle, the positron. In a hydrogen atom, the Bohr radius and Rydberg energy are 

5.29 x 10711 m and 13.6 eV. You may take the reduced mass of a hydrogen atom to be the same as the mass of an 
electron. 


Topic 7 — Relativistic hydrogen atom 


Exercise 7.1 The Coulomb model of a hydrogen atom treats the electron and proton as point charges that 
interact only via the Coulomb force. List the additional effects that occur in real hydrogen atoms and give rise to 
fine and hyperfine structure in the hydrogen atom spectrum. 


Exercise 7.2 (a) In the absence of the spin-orbit interaction, the hydrogen atom has 8 degenerate quantum 
states with n = 2. Specify these states in terms of their n, l, m and ms quantum numbers. 


(b) When the spin-orbit interaction is included, the n = 2 states are classified by the quantum numbers n, l, j 
and mj. Specify the 8 hydrogen atom states with n = 2 in terms of their n, l, j and m; quantum numbers. Which 
of these states have the same energy according to relativistic quantum mechanics (which includes spin-orbit 
effects)? 


Exercise 7.3 In the first-order perturbation theory of the relativistic hydrogen atom, every unperturbed energy 
level is affected either by the spin-orbit interaction or by the Darwin term, but not by both. Explain why this is so. 


Exercise 7.4 (a) Express L - S in terms of i. g” and T = (L + s)’. 


(b) Consider two quantum states, |A) and |B), in a hydrogen atom, both labelled by the same values of n and J, 
with / > 0, but with j = l + 5 and j = l — 2, respectively. Use your answer to part (a) to show that the first-order 
energy corrections for these states, due to the spin—orbit interaction, are in the ratio 


EO ıl 
pO THI 


Topic 8 — Many-electron atoms 


Exercise 8.1 Explain why, in the central-field approximation for a many-electron atom, energy eigenfunctions 
can be labelled by the quantum numbers n, l and m, just as in a hydrogen atom. 


In the Coulomb model of a hydrogen atom, the energy levels depend on the principal quantum number, n, but are 
independent of l. Why does this not apply in many-electron atoms? For a given value of n, does the energy 
increase or decrease with |? Give physical reasons for your answer. 


Exercise 8.2 (a) List all the atomic shells for n = 1, 2, 3 and 4, and give their capacities. 


(b) Given that the energy ordering of atomic shells is: 1s, 2s, 2p, 3s, 3p, 4s, 3d, ..., write down the ground-state 
configuration of a silicon atom (Z = 14) and an iron atom (Z = 26). 


Topic 9 — Atomic terms and levels 


Exercise 9.1 An atom has a pair of valence electrons in open shells with (l = 1, sı = 5) and (l2 = 3, sg = 5). 
Specify the atomic terms that arise from this configuration in terms of the values of their L and S quantum 
numbers, and give an appropriate spectroscopic symbol for each term. 


Exercise 9.2 The excited states of a helium atom have one electron in the lowest energy atomic orbital and the 
other electron in a higher energy atomic orbital. In these excited states, triplet states have lower energies than 
singlet states. Explain why is this so, using an argument based on first-order perturbation theory. Your argument 
need not contain any equations, but should emphasize the role of identical particles and electron—electron 
repulsion. 


Topic 10 — Molecules 


Exercise 10.1 A beryllium atom, Be, has atomic number Z = 4. Write down the ground-state electronic 
configuration of a diatomic beryllium molecule, Beg, and determine its formal bond order. You may assume that 
molecular orbitals of Beg are in order of increasing energy: log, loy, 20g, 20u, ..., where lo, and 20g are 
bonding orbitals and lo, and 20, are antibonding orbitals. Would you expect a Bez molecule to have a high 
dissociation energy? 


Exercise 10.2 (a) What is the physical basis of the Born—Oppenheimer approximation? How does the 
assumption of the Born—Oppenheimer approximation allow the time-independent Schrödinger equation for a 
molecule to be simplified? 


(b) What is an energy curve for an electronic state of a diatomic molecule and how does the assumption of the 
Born—Oppenheimer approximation allow such a curve to be constructed? 


Exercise 10.3 (a) Consider a hydrogen molecule in the Born—Oppenheimer approximation. Suppose (for the 
purposes of this question) that the ground-state electronic energy eigenvalue is given by 


2 4 
ë la 2 
EaR) = — (75-5), 
4reo \5R R 
where a is a constant length and R is the proton-proton separation. Write down an expression for the total static 


energy of the molecule in this electronic ground state. 


(b) Use your answer to part (a) to determine the equilibrium separation of the protons and the spectroscopic 
dissociation energy for a hydrogen molecule in its electronic ground state. 


Exercise 10.4 In the LCAO approximation, the ground-state energy eigenfunction for a hydrogen molecule ion 
has the form 


A B 
y(r) =C ( is(r) ale oh (r)) , 
where o and oR are real-valued atomic orbitals centred on the two protons, A and B, and C is a normalization 
constant. 
. . . . . . A B 
(a) Derive an expression for C in terms of the interatomic overlap integral (¢;\|¢7,)- 


(b) Explain why ~;(r) corresponds to a bonding molecular orbital, while a function of the form 


yal) = D (lr) — oR (x) , 


where D is a suitable constant, does not. 


Topic 11 — Solids 


Exercise 11.1 State Bloch’s theorem for the state of an electron in an infinite crystalline solid. Show that it 
implies that the electron probability density is the same at all lattice points in the crystal. 


Exercise 11.2 Consider the function 


u(r) = 5y o(r — R,)eik(Bi—) 


where each R; is a lattice vector, the sum is over the lattice vectors in an infinite crystal and k is a real constant 
vector. 


(a) Show that u(r) has the periodicity of the lattice. 
(b) Hence show that the function 


Ye) = D or — Re, 


where each R; is a lattice vector and the sum is over the lattice vectors in an infinite crystal, satisfies Bloch’s 
theorem. 


Exercise 11.3 A donor atom in its lowest energy configuration in germanium has an extra bound electron 
compared to a germanium atom. How much energy must be supplied to the donor atom for the extra electron to 
become unbound? You may treat the donor atom plus its additional electron as a hydrogen-like atom with an 
effective atomic number Z = 1 and assume that the effective mass of an electron in the conduction band of 
germanium is mz = 0.12me¢. The relative permittivity of germanium is £ = 15.8 (making the Coulomb force 
between charges in silicon smaller than that in empty space by a factor 1/15.8). 


Exercise 11.4 Explain why the electrical conductivity of a typical metal decreases with increasing temperature, 
while the electrical conductivity of a typical pure semiconductor increases with increasing temperature. 


Topic 12 — Interaction of light and matter 
Exercise 12.1 (a) Which types of radiative transition have transition probabilities that can be predicted 
directly using first-order time-dependent perturbation theory? 


(b) If stationary states are timeless and unchanging in quantum mechanics, how is it possible for a transition to 
occur spontaneously between two stationary states in an isolated atom? 


Exercise 12.2 For monochromatic light polarized in the z-direction, the electric dipole approximation introduces 
a time-dependent perturbation 


VG@)={@2= FU) r cose, 


which is independent of the azimuthal angle, œ. Show that, according to the results of first-order time-dependent 
perturbation theory, such light will cause no radiative transitions between two hydrogen atom states that have 
different values of the magnetic quantum number, m. 


Exercise 12.3 For monochromatic light polarized in the x-direction, the electric dipole approximation introduces 
a time-dependent perturbation 


V(t) = f(t)x = f (t)r sin 8 cos ¢, 


Show that, according to the results of first-order time-dependent perturbation theory, such light will cause no 
radiative transitions between two hydrogen atom states that have the same value of the magnetic quantum number, 
m. 


Exercise 12.4 (a) Define the three different types of radiative transition. 
(b) The following equation applies to atoms in equilibrium with radiation: 
Bizu(wa)Nı = Bau(wz)N2 + A2102 


where Bj, Bo; and Av; are constants, N; and No are the numbers of atoms in quantum states 1 and 2, of energies 
E and E> (with E2 > E1), u(w) is the spectral energy density function and w21 = (E2 — E1)/h. 


Correlate the three terms in this equation with the three types of radiative transition defined in part (a). 
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Solutions 
Topic 1 — Spherical harmonics 


Solution 1.1 (a) We have 
P 2 
Tro= -r| o (noZ) "E. =| Acos 


sin 6 06 00 sin? 0 0¢? 
1 d dcos 0 
—_ apl tfa 
=e ode (sno dd ) 
i d 
_ Ape a 7.2 
SA gag 
= Ah? q x 2 sin 0 cos 0 
sin 0 
= 2h? A cos 8, 


so f (0) is an eigenfunction of L’ with eigenvalue 2h. 
Also, 


L.Acos@ = ih A cos 8 


= 0 
= 0 x Acos8, 


so f(@) is an eigenfunction of L, with eigenvalue 0. 
(b) We require that 


2T fT 
rs f f |A|? cos? @ sin 0 dé dọ 
0 0 


= arak f cos? 6 sin 0 d8 
0 


2 
= 27| Al? x 3 


4 
= =n Al’, 
3 


so we can choose A = 1/3/47. 


In general, the eigenvalues of T” and L, are I(l + 1)h? and mñ respectively, so (4/3/47) cos 0 is a normalized 
spherical harmonic, Yım, with quantum numbers l = 1 and m = 0. 


Comments (1) There is no problem with having an eigenvalue equal to zero. Whenever Ô f(x) = 0, we can say 
that f(x) is an eigenfunction of O with eigenvalue 0. 


(2) Normalization of an angular function f (6, ġ) requires that 


2m T 
f 1 |f(0,6)|? sinf dð dọ = 1. 
0 0 


Note that the integrand contains a factor sin 0 and the limits are: 0 < 0 < m and 0 < ¢ < 2r. 


Solution 1.2 (a) Using Dirac notation, the angular part of the wave function is 
lang) = $|1,1) + 511,0), 
so 
(anglang) = [2(1, 1| + 4(1,0|| [2|1, 1) + 2|1,0)| 
= (1, 11,1) + $2 (1, 1/1, 0) + 32(1, 0/1, 1) + # (1, 0|1, 0). 
11 


The spherical harmonics are orthonormal so (1, 1/1, 1) = (1, 0]1,0) = 1 and (1, 1]1,0) = (1,0|1,1) = 0. 
Hence, 
(anglang) = 95 + 55 = 1, 
as required. 
(b) Both parts of the wave function contain spherical harmonics with | = 1. The wave function is therefore an 


a2 
eigenfunction of L with eigenvalue 1(1 + 1)h? = 2h?. The probability of getting the value 2h? in this state is 
equal to 1. 


(c) The angular part of the wave function is normalized so we can use the coefficient rule to read off the 
probabilities of the possible values of L}. The probability of getting the value A is found by taking the square of 
the modulus of the coefficient of |1, 1), which gives |3/5|? = 9/25 = 0.36. 


Comments (1) It is convenient to use the Dirac notation |l, m) for Yi;,(0, ¢), as it avoids the need to write down 
explicit integrals. Of course, 


20 T 
(astaltasma) = | f Yim (0.0) Yima (0,0) sind 40.40. 
(0) (0) 


(2) We are told that the whole wave function is normalized and we have shown that the angular part is normalized, 
so it follows that the radial part is also normalized: 


T |f)? r? dr =1. 
0 


Because of this, there is no need to worry about the radial part of the wave function when using the coefficient 
rule. Remember, the coefficient rule comes from the overlap rule. If the overlap rule were used directly to 
calculate the probability, the radial function would integrate out to give a factor of 1. 


Topic 2 — Total angular momentum 


Solution 2.1 Denoting Yj,,(0, Q) by |l, m}, we have 


|) = 31,011) — 2 1.) Y) 
|B) = /2 11,0) 1) + 411,14). 


The spherical harmonics are orthonormal, so we have 


27 T 
(l1, millo, m2) = f | Yi m: (9, ®) Vine (0, $) sin 6 dO do = ôi, lı Ôm2,mı- 
o Jo 


Hence, 
(Ala) = (4/8 Lot- y3 ae NLY - 2 1,14) 
= 4(1,0/1,0)(f | t) — (1, 011,1) 14) 
— PULOU IY + 31, 11,1) |) 


=4(1x1)- 20x 0)- Box 0) +2 (1 1) 
1 


(BIB) = (V3 Lor l+ VE e N (VRLO D + VFI, 01) 
= 2(1, 0/1, 0) (t | t) + ELLIN 
+ 2 (1, 1/1, 0) | 4) + $1, U1, 1) |) 


=2(1x1)+20x0)+ 20x 0)+40x1) 
1 


12 


So |A) and |B) are both normalized. Also 
(BIA) = (9/2 (oye | + yE ald 1) (4/4 11,091 4) — 9/2 11,091 )) 
= (1, 0/1,0)(t | t) — 3(1, 0/1, 1) (4 |) 
+ (1, 11,044 | t) — (0, 12,1) 44 | 4) 


(1 x 1) -— 2(0 x 0) + 4(0 x 0) - Bx 1) 


AS 


= 0. 
Hence |A) and |B) are orthogonal. 


Comment The space and spin parts of the total wave function act like they don’t see one another. So, when we 
form the inner product of the total wave function and multiply out terms, the spin parts come together to produce 
their own inner products and the space parts do the same. The final answer is therefore a sum of terms, each of 
which is a product of two inner products — one involving spin states and one involving spatial states. 


Solution 2.2 With 


= (1/2 Mol DF J3¥ial p) 
we have 
J.|B) = (L. + 8.)|B) 

= (L. +8.) (1/2 Mol T) + VE Mil \)) 
= /3(E.%10)1 D + 3 (EY) + y4 Yio(8:1 0) + y4 Yia (8.1) 
= 05/2 iol t) tAE Mal 4) + day’? Mol t) — tayi Nal) 
= $n(/2%.01 4) + /3Mial 4) 
= 3ħ|B), 

so |B) is an eigenvector of J, with eigenvalue h/2. 


Comment In an expression like (L z + S.)¥io| +), we multiply out terms and remember that L z acts only on the 


spatial part (Y1 ,o in this case), while S; acts only on the spin part (| t) in this case). (To keep the notation compact 
we have omitted the angular arguments of the spherical harmonics, but we must not forget that these functions 
depend on the spatial coordinates 6 and ¢, and so are affected by the operator Lz.) 


Solution 2.3 (a) For |l, j,m ;) = |1, 5, $), we have 


~~ 2 <2 
L- S| aa) 5 (J -L -SLAI 
1 
= (a3 +) -1+1 -3G +1) F711, 3,5) 
= —h*|1, 5, 4). 


For |l, j, mj) = |1, 3,5), we have 


Get) 10+) 304 1)r las) 


z ia, $4) 
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Hence |1, 5, 5) and |1, 3, $) are eigenvectors of L - S with eigenvalues —h? and +h?/2, respectively. 


(b) With spin-orbit coupling, |1, 5,5) and |1, 3, $) 
configuration in which the electron has / = 1. The spin-orbit interaction is proportional to L- S, so the fact that 


these two eigenvectors have different eigenvalues of L - S leads to spin-orbit splitting in l = 1 configurations. 


are atomic terms that develop from an unperturbed 


Comments (1) By definition of the scalar product, 
L- S= L + L,Sy+ ea 
but this is a very inconvenient expression for calculations because Ly, Ly and L, cannot all have definite 


non-zero values in the same state (and similarly for S,, Sy and S,). Nevertheless, the particular combination 


LS, + LS, + LS, does have a definite value in any state with definite values of J?, L? and S?. We can find 
this value by using the identity given in the question. 


(2) Note that L - S is the same as § - L because each of the spatial operators Lr, Ly and Le commutes with each of 
the spin operators Sz, S, and S.. 


Topic 3 — Hydrogen atom 


Solution 3.1 In the Coulomb model of a hydrogen atom, the energy levels are 


where n = 1, 2, 3, ..., so the energy level — Ær /16 has n = 4. The possible values of l are therefore 0, 1, 2 and 3. 
For / = 0, the only allowed value of m is 0. 

For / = 1, the allowed values of m are —1, 0 and 1. 

For | = 2, the allowed values of m are —2, —1, 0, 1 and 2. 

For l = 3, the allowed values of m are —3, —2, —1, 0, 1, 2 and 3. 


This gives a total OF 1+3-+5+7 = 16 combinations of l and m for n = 4. For each of these combinations we 
can have m, = +5, so there are 16 x 2 = 32 degenerate states with energy — Eg /16. 


+4, 


Comment In the Coulomb model of a hydrogen atom, 2n? degenerate states are associated with the quantum 
number n. For n = 4, this gives 2 x 4? = 32 states. This degree of degeneracy is lowered in the fine structure of a 
hydrogen atom (where energy levels depend on the n, / and j quantum numbers) and in many-electron atoms, 
even in the central-field approximation. 


Solution 3.2 Using the method of separation of variables, the time-independent Schrödinger equation for the 
Coulomb model of a hydrogen atom can be split into separate equations for the r-, 0- and ¢-dependencies, 
allowing us to write 


Wrtrall 6, p) — Rni(r)Oim(0)e""?. 


The magnetic quantum number m must be an integer to give 
eim(o+2r) = eme, 


This is needed to ensure that the energy eigenfunction is single-valued in the sense that 
Wnim(T, 0, Q a 27) _ Wnim (r, 0, $). 


The function O©;m (0) diverges as 0 — 0 and 0 — ~ unless | is a positive integer greater than or equal to |m|. The 
function Rw (r) diverges as r — oo unless n is an integer greater than l. These quantum numbers take integer 
values because it is assumed that the energy eigenfunction does not diverge anywhere. 


Comment The requirement that the energy eigenfunction be single-valued is not entirely obvious. Certainly, 
ly(r, 0, @)|? must be single-valued because it represents a physically-measurable quantity — the probability 
density of finding a particle at a given place. The single-valued nature of (r, 0, @) is less clear-cut because an 
eigenfunction is not a measurable quantity. Nevertheless, once an overall constant phase factor has been chosen, 
quantum mechanics assumes that eigenfunctions and wave functions have definite values at each point in space, 
and this assumption gives the quantization condition for the magnetic quantum number m. 
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Solution 3.3 Using the given form of radial function, the radial probability density is given by 
Prins (Crr-te-r/na0)’ = C?r eT? /nao, 


This is a maximum when 


0= < (Ce) 


2n 
= C2 | Inr?! — 2r e7 2r/nao 
nag 


= r = 
E 2022" 1 ae e 2r (m. 
nag 


which gives r = nap. 


For l = 99 and n = 100, the typical size of a hydrogen atom is 1007a9 = 5.29 x 10~7 m (comparable with the 
wavelength of visible light). The energy of an atom in this state is —13.6 eV /100? = —1.36 meV, so the energy 
needed to ionize the atom is 1.36 meV. 


Comments (1) Note that the radial probability density does NOT contain a factor of 47. This can be traced 
back to the fact that Ynim(r, 9,6) = Rni(r)Yim(@,¢), where the spherical harmonic Yim (0, Q) satisfies the 
normalization condition 


27 T 
f | \Yim(0,)|?sin@d0dgé=1 (NOT 47). 
0 0 


(2) Hydrogen atoms in states with very high values of n (Rydberg states) have been prepared experimentally, but 
they are very easily ionized. 


Solution 3.4 The expectation value of r in the given 2s state is given by 


(r) = I Ry 0(r) r Roo(r) r? dr 


1 j 2 r ? eee he 3q 
= aap ff @-a) leo) te 


1 love) ; 4 4 5 
= ox | Gi Ta + =) e7™"/% dr, 
(2a0)? Jo ao aĝ 


Using the standard integral given in the question and noting that 3! = 6, 4! = 24 and 5! = 120, we obtain 


1 24a 120aĵ 
= — |4 x 6ag — 4 x — 2| = bao. 
(o) = gag [A x 6af — a x 08 + E — ba 


Comment If we know the radial function R,ı(r), the expectation value of any function of r can be found 
without doing any angular integration. In any state described by the radial function R,,; (7), the expectation value 
of f(r) is given by 


(f(r) = f "pla 9? dr. 


There is no factor of 47 in this expression for reasons similar to those given in the preceding question. 
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Solution 3.5 (a) We have 


il 20 T lee) . 1 l 
(=) =| / f P21.0(r,0, >) -z3 V21o(r, 0, ġ) r? sin 6 dr dé do 
2m 2 
= (5 ya A f [= = (reni cos J r? sinð dr dé dọ 
~ Ot 2a9 az 
1 3 27 j 
= —r/fao 
ie ~ (5) a ao f COs *osinoae |” re dr. 


Using the standard integrals given in the question, we obtain 


1 1/1 i oe te? 1 
— = — — T = an = =. 
r3 4n \2ao) a? 3 24a3 


(b) The spin-orbit interaction is proportional to 1/r?, so finding the expectation value of 1/r° in various states 
of a hydrogen atom is needed when calculating the first-order perturbation approximation for the spin-orbit 
splitting of hydrogen atom energy levels. 


Comments (1) In this question the complete energy eigenfunction Ynim(r, 0, ¢) is given, but there is no obvious 
way of splitting this into a radial part and an angular part, each with appropriate normalization. We therefore resort 
to the sandwich rule for expectation values, using the complete energy eigenfunction and integrating over 0 and ¢ 
as well as r. 


(2) As in all volume integrals in spherical coordinates, we must use the volume element dV = r? sin 6 dr dé dd. 


Solution 3.6 (a) The expectation value of r is given by 


27 
»= [ I r|Ral(r)Yim(0, 0) |? r? sin 6 dr d8 de 
oo Qn T 
7 r8|Ru(r)Pedr | i lYim(0, b)|? sin 8 dé do 
0 0 0 


= f r?|Rni(r)|? dr, 
0 


because the spherical harmonics Ym are normalized. Similarly, 
OO 
(2) = f Ratar, 
0 


so (r) and Ar = 4/ (r?) — (r)? are independent of m. 
(b) The probability density at any given point is 
lVnim(r, 0, p)? = |Rr (r) Yim (4, 6)/? 
= |R (r)? |Om (8)? le? |? 
= [Ru (r)|’lOim (0), 
which is independent of ¢. 


(c) The radial probability density at any given radial coordinate is given by r?|R,1(7r)|*, which is independent of 
m, 0 and ¢ (by definition, the radial probability density is a function only of r). 


Comment Part (a) could be answered more quickly using the general formula 


= [FORO 
0 


but the above method shows explicitly why the dependence on m disappears when the angular integrations are 
carried out. 
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Solution 3.7 We have 


Ble 
gle 
N 
bo 
Sja 
| 
4 
Pa, 
2 
L=] 
pO 
| 
3| = 
om 
Ç| e 
mes 
S13 
© 
| 
af 
bisa 
a 
(æ; 
eee 


a2 a2 
The operator L” involves partial derivatives with respect to 6 and œ, and so L e~"/% = 0. Hence 


pre tis — — h2 V2e7"/20 


and so 


PL R 2ao 
PY, = = (= = i) W1,0,0- 


2 
a 
We therefore have 


2T T lee) 
= / I f Vto.0(7, 9,9) DP W1,0,0(7,0, ġ) r° sin 6 dr dé do 


27 
A A s WPi o,o(r, 8 9) a = (#2 — i) W1,0,0(7, 9, Q) r? sin 6 dr dé do 


27 
_ Fi [s e7r/ao 2 ne- 1)e =r/% r2 sin ð dr dd dd 
nas að id 


R2 
= — e (2aor — r-\e =2r/a qr 
0 


= 
ag 


Comment The expectation value of any quantity has the units of that quantity. Recalling the de Broglie relation 
(AaB = h/p), and noting that ao is a length, it is easy to see that our answer has the correct dimensions for the 


expectation value of the square of momentum. 


Topic 4 — Perturbation theory 


Solution 4.1 Choose the unperturbed Hamiltonian to be that for a particle of mass m in the corresponding 
one-dimensional infinite square well, with 


0 for —L/2 < x < L/2 
ee) elsewhere. 


In the ground state, the unperturbed energy eigenfunction is 


with eigenvalue gO = n?ħ?/2mL?. The perturbation is then 
A Alz] for —L/2 < x < L/2 
ôH = 


0 elsewhere. 


So first-order perturbation theory gives a first-order correction 
1 0)) fF] (0 
EP = (WP i). 


Because the perturbation involves |x|, the integral we need to evaluate has an even integrand, and we can write 


We change the variable of integration to u = (7a/L). Then, x = (L/7)u and dx = (L/7) du. Also x = 0 and 
x = L/2 correspond to u = 0 and u = 7/2, so 


2 pr/2 
ED = 4A (=) f u cos? udu 
L \r 0 


-Afi  W—4 
OL \er 16 
AL 


= qa E 4). 


The first-order approximation to the ground-state energy is 
nh? AL 


By ~ ——~ + —> 
1 mL? | 4x2 


(x? — 4). 


In the first excited state, the unperturbed energy eigenfunction is 


yw) n= = sin (=) 


with eigenvalue EW) = 4n*h? /2mL?. 


So first-order perturbation theory gives a first-order correction 


EP = YP oH |W?) 


The first-order approximation to the energy of the first excited state is 


Ar?h2 AL 


IAL? j 4 ` 


9 = 


Comment For À > 0, the energies of both states are increased by the perturbation, but the increase is smaller for 
the ground state than for the first-excited state. This is because the probability density in the ground state is largest 
around x = 0, where the perturbation is equal to zero. By contrast, the probability density has maxima away from 
x = 0 in the first-excited state, where the perturbation is larger. 
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Solution 4.2 Choose the unperturbed Hamiltonian to be that for a particle of mass m in the corresponding 
one-dimensional infinite square well, with 


0 for —L/2 < x < L/2 
VO(xr) = 
ee) elsewhere. 


In the ground state, the unperturbed energy eigenfunction is 


2 
Pa) =y 5008 (>) 
with eigenvalue gO = n? h? /2mL?. The perturbation is then 
de A for —w/2 < x < w/2 
ôH = 
0 elsewhere. 
The first-order correction to the ground-state energy is 


EO = (yP iy) 


= z [doo (=) dz. 


Change the variable of integration to u = mg/L. Then dx = (L/7)du. Also, x = +w/2 corresponds to 
u = +rw/2L. So 


QL mw /2L 
EW = == | cos? udu 
La —7w/2L 


7 2) E 1 _ 


=+ A sin(2u) 


qm |2 


—rw/2L 


_ 4X eee (=) 
el ao NE 


=a) P+ Esin (Z8)]. 


The ground-state energy of the system is therefore estimated to be 


‘a w 1 TW 
Ex pa nme (=) 
1 Ome [Z+ sn L | 


Comment The correction tends to 0 as w — 0 and tends to À as w —> L, as might be expected. 


Solution 4.3 We use perturbation theory with the potential energy function split into an unperturbed part 1740) 
and a perturbation V\) where 


e2 e7"/à 


2 
1 ps 
al: 2 ad A= 
(r) ATEQ Tr z 4neo r 


Perturbation theory then gives the first-order energy correction to be 


E® = (4 0,0|6Ei|¥1,0,0) 
2T T ee) 
= | | f Wio,0(7, 0, 6) ôH Y1 ,o,0(r, 0, 6) r° sin 8 dr dd do 
o Jo Jo 
2 oo —r/X 
T = (- £ ) ar | e7r/a & e77/% p? dr 
Tag Are 0 r 


2 


e 0O 
5 a e- (2/a0+1/A)r p» dy. 
TEQ AG 0 
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Using the standard integral given in the question, we obtain 


2 
Bo) œ 1 
to TEQAg 2 1 2 
(a +4) 
e A 


TEQ ag (ao + 2A)? ` 


Comment Without the perturbation, the ground-state energy in the Coulomb model can be expressed as 
2 


EY =-= -— 
1 i 8T E040 i 
so the ratio of the first-order correction to the unperturbed ground-state energy is 
E” B e? N si 8TEQAg 8A? 
EO E TEQ ao (ao + 2X)? e2 7 (ao + 22)?’ 


and for À < ap, this is approximately equal to 8A?/ Gn So, if the ground-state energy agrees with the Coulomb 
model to one part in 1000, À cannot be much larger than y aĝ /8000 ~ 0.01ao. This provides a check on the 
validity of Coulomb’s law at short distances. 


Solution 4.4 A hydrogen atom is a two-particle system characterized by a reduced mass u. We can use the 


general formula for the relativistic perturbation to the kinetic energy (given in the question) with m replaced by p. 
The first-order correction is then 


i pS E X 
Ey” = (trol 5Hl1.00) = - a (P Viool 00), 


where u is the reduced mass for a hydrogen atom. Hence 


2r 
DROL L COS Pe) man 
2r 2 
leh LL GEJAN GEJA) ema 
mag \ 8u3c ag Xt 
— 1 4ao dag e7 27/40 r2 
E o a aaa 


3 
0 
ht love) 
=~ ay f (r? — 4agr + dap) e72r/% dr, 


BP = raag EOD -te (P) +408 (2) 


So 


where Er = h?/ 20a? is the Rydberg energy. 


Comments (1) Because Ep and uc? have the dimensions of energy, our final answer for the energy correction 
has the correct units of energy. 


(2) The trick of using the Hermitian character of p° to write (|p ly) = (P Yp y) is a useful one because it 
avoids having to work out the fourth derivative. In fact, the calculation of the fourth derivative is hazardous in this 
case because there is a singularity at the origin. This can be handled using an advanced mathematical technique 
beyond the scope of this course (the so-called Dirac delta function) but it is good to avoid such complications. 
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Solution 4.5 We have 
Hy|Un,a) = Ey las): 


SO 


(Yn, slHlYn,a) = En, A Unan, a) = Ena; 
because the eigenvectors are normalized. 
Differentiating both sides gives 


“Fe = (e 


iaa, 


h 


Yn a)+ (una lea 


+ [tn a) + (Wn, A 


Using the fact that f à is Hermitian, and that y,,,, is an eigenfunction of H à with the real eigenvalue E,,,), we have 


(2h > Un a= En al Pea 3) 
and 
(vna Ba A) = (Pona G54) = Ena (tna A) 


Adding these two equation then gives 


fi, |e = Sii (ee 


(“enol 


Un a) + (un, A 


Yn a) + (vn, A 


a) 


En, È lna 


=0; 


nr) 


because (Yn a|n, à) has the constant value of 1 irrespective of A. Combining this with our previous expression for 
OE;,,,/OA, we see that only the second term on the right-hand side of this equation survives and so we obtain 
EnA 
Or 


as required. 


dH) 


= (Ynal S 


Yn as 


Comments (1) Because the energy eigenfunctions depend on the parameter A, the derivative of the expectation 
value in principle involves derivatives of the eigenfunctions. The key point about the Hellmann—Feynman theorem 
is that these derivatives cancel out, and we are left only with the expectation value of the derivative of the 
Hamiltonian operator with respect to A. 


(2) When using the Hellmann—Feynman theorem we generally omit the subscript À in our notation for the 
Hamiltonian operator and the energy eigenvalues and eigenfunctions. It is included in the above proof to make 
the dependence on A explicit, but we have no need to include it when applying the theorem — we just take the 
parameter À to have its physically appropriate value. 


Exercise 4.6 (a) From the given Hamiltonian operator and energy eigenvalue, we have 


OH, Ces iai te  m 
—_ >= = d == = h = = 
ag 22 o a oe 20” 
so the Hellmann—Feynman theorem gives 
E, 
al = (Un| 5x? ln) = (52°). 
Multiplying both sides by C then gives 
5 En = (402?) = (Epot); 


as required. 
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(b) From the given Hamiltonian operator and energy eigenvalue, we have 


oH 2e ðEn 4 pi En 
— = — ani = — —_ — = — 
ðe Areor de (4re0)? 2h? n? e 
so the Hellmann—Feynman theorem gives 
En 2e 
7 = (Ynim| = Gmeor (Ym): 
Multiplying both sides by e/2 then gives 
e? 
2En = (YPnim| = Taar Vm) = (Epot) 


as required. 


Note: Since F, is negative, this implies that (Epot) is twice as negative as the energy eigenvalue. This is 
compensated by the fact that (Fin) has the positive value —E,,, so (Ekin) + (Epot) = —En + 2Ey = En. 


Comments (1) The Hellmann—Feynman theorem is not a formal part of the course, so you need not remember 
it. As you have seen here, the theorem is useful if differentiating the Hamiltonian operator with respect to a 
parameter gives another operator whose expectation value in an energy eigenstate is of interest. We need to have 
an explicit formula for the energy eigenvalues. 


(2) Remembering that the total energy is Exin + Epot, the results derived in this question can be re-expressed as 
follows. In any energy eigenstate of a one-dimensional harmonic oscillator, (Ein) = (Epot), and in any energy 
eigenstate of a hydrogen atom in the Coulomb model, (Ekin) = — (Epo). The latter result is possible because the 
expectation value of the potential energy is negative in this case. 


Topic 5 — Variational method 


Solution 5.1 (a) We have 


ks — are 2 
x 
so 
q? —aa? /2 2,.2\ ,—ax~/2 
ae =(-a+a*x*)e 
Hence, 
h2 lee) love) 
(h T Tai H) TA (0 D ee” dr — a? a reo?” ar) 
(dil Pt) I e702? dz 
h2 Ro, L ae dr 
om m |e 2 
f e dz 
oi an l 
7 Im” 2m 2a 
h2 
= fon 
Also, 
L 2 
e "dx 
CAC Í 
(pilo) 


29 2 
f e ° dz 
-00 
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Change the variable of integration to u = yaz, then 


Jal 
=u? du 
[ i Vo [ME 2 


(AVe) _ _y, 2 et ar du. 


(ile) 


So 


(afle) h? Ve [~VE 
Re = —a- -2 du, 
Od) m° a _ 


as required. 


(b) When a is positive, but very small, we can approximate the integral by 


val, 
f e™ du ~ e x VaL = vaL. 
0 


So, in the limit of small positive a, we have 


h2 
4m” rt avaL, 


and this is negative when 


VoL 
X fa Z 
VT’ 
i.e. when 


E, = — 


0<a< W 
ħin ’ 
which can always be achieved. The exact ground-state energy must be smaller than any variational estimate, so the 
exact ground-state energy is always negative (below the top of the well) no matter how small Vo and L may be. 


Comment We have shown that any finite square well in one dimension (no matter how narrow or shallow) has at 
least one bound state. This result was quoted in Book 1, but not proved there. However, this result turns out to be 
specific to one-dimensional wells. It is possible to have finite wells in two or three dimensions that are so narrow 
and/or shallow that they have no bound states. 


Solution 5.2 Using the identity given in the question, we have 


h2 2 h2 1 d 2q —r/a 
Sam V8) = ~ dr ¢ ar® ) 


h? 1 d r? —r/a 
= —— — | —e 
2mr2?dr \a 


Fl w e t/a 
2mr? \a a? 


Thus, 


We therefore have 


27 1 1 
(|b) = il [ I om (=- =) + gmbr”? e '/*r? sino dr dO dd 


h? 1 a4 2 
=s f =. sa + z0r | e72"/a dr. 
0 


Using the standard integral given in the question, 


ion =e (at G) — aa G + (3) 


Ra Ra 3 


We also have 


27 T love) 
(Qdo) = i. i f e 27/4 7? sin A dr dO do 
0 0 0 


= tn | re 2"/4 dr 
0 


= Ar 2! (5 
2 


= 7a", 
so 


f h2 
E= (GIAO) _ P muda? 


(Plot) 2ma? 


The minimum value is obtained by setting 


dEr h? 
j= da = E? 3 + 3mwĝa, 
SO 


and the estimated ground-state energy eigenvalue is 


ie h 


Ermin = — 
een 2m h 2 V/3mu 


Comment This estimate is about 15% greater than the exact ground-state energy (see next exercise). The 
variational method always produces an estimate that is greater than, or equal to, the ground-state energy. 


Solution 5.3 We have 


— B 
-5 VAr) = -5 
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We therefore have 


~ ae fE fS jo, | K r? 1 2 19,2 
(d|H| dé) = f f f ec | (s — =) + ymar?! eo” /20" 5? sin 9 dr d8 do 
0 Jo Jo 


2ma? 


= ar f por ~a nya e EG dr. 


Using the standard integrals given in the question, 


Bh? YT 3 B 3V7 5 
Ima? 4° mal 8 “ 


ON 


lallo) = 4r ( Ta) 


1 
2 
R? h? 1 
= 127./7 (E = Tema + ree) 


_ 3ayT 6 r ) | 


—a+mwa 
4 m o 


We also have 


27 T love) 
(lb) = I f f e” /” 72 sin dr dd d 
0 0 0 


= an | re /@ dr 
0 

= 4r VT 
4 

= ny Ta’, 


so 
(afle) T p? 2 i) 
k= — = -| — + mwĝa’ }. 
© (o) 4 Xa? © 
The minimum value is obtained by setting 
dz, 3 f 2 3 
SE Sf Sl D 
da 4 ( tr a 
so 
ħ 
go 
MWO 


and the estimated ground-state energy eigenvalue is 


3 (K mw h 3 
Ermin = 7 =— + map) = zwo: 


Comment This estimate turns out to be exact because the trial function e~” /” happens to include the exact 
ground-state energy eigenfunction (which has a? = h/mw0). The estimate made in Exercise 5.2 is slightly greater 
than this, because a trial function of the form e~"/* is less accurate for this system. 


Topic 6 — Hydrogen-like systems 


Solution 6.1 For helium-4, we have m = 7360m, so the reduced mass is 
207 x 7360 

= Me 
207 + 7360 


We are told that the reduced mass of hydrogen can be taken to be me, so 4/44 = 201.3. For helium, Z = 2, so the 
scaled Rydberg energy is 


u = 201.3me. 


Egad _ 72 É Fp = 92 x 201.3 x 13.6eV = 1.10 x 104eV 
HH 


and the scaled Bohr radius is 


1 
gaed. = ZI” = 5.29 x 1071 m/(201.3 x 2) = 1.31 x 107! m. 


Comment This calculation assumes that the muon is mainly outside the helium-4 nucleus, which is justified 
because the nuclear radius of helium-4 is of order 1.7 x 10715 m 


Solution 6.2 The reduced mass of positronium is 
B m2 _ Me 
i Metm 2 
We are told that the reduced mass of hydrogen can be taken to be Me, so u/ {uy = 1/2. The scaled Bohr radius is 


PH a0 = 2ao = 2 x 5.29 x 107! m = 1.06 x 107! m, 
u 


and the scaled Rydberg energy is 
i 
a É fp = 5 X 13.6eV = 6.8eV, 


HH 
so the ground-state energy is 


psealed 
E=- 2 = -EI — _6 Bev. 


Comment The usual formula for the scaled Rydberg energy in an atom involves the atomic number Z, which 
appears because the charge of the nucleus is Ze, where e is the magnitude of the charge of an electron. In 
positronium, the positron plays the role of the nucleus in an atom, so in this case Z = 1. 


Topic 7 — Relativistic hydrogen atom 


Solution 7.1 The additional effects include: 


e A relativistic correction to the kinetic energy 
e the spin-orbit interaction (for l Æ 0) 

e the Darwin term (for l = 0) 

e the Lamb shift 

e the finite size of the proton 


e the magnetic interaction between the magnetic dipoles of the electron and the proton. 


Comment The first three effects give rise to fine structure (corrections of order a”, where a ~ 1/137 is the fine 
structure constant) while the last three effects give rise to hyperfine structure. 
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Solution 7.2 (a) Inthe absence of a spin-orbit interaction the states are: 
n, l, M, Ms = 2,0,0, +5) |n, l,m, ms) = |2,0,0, -3) 


n, l, M, Ms = 2,1,1, +5) |n, l,m, Mms) = |2,1,1,-3) 


) 
) 

tLe, Ms) = 2,1,0; +4) |n, l, m, Mms) = |2, 1,0, —4) 
) 


n,l,m, Ms) = 2,1, —1,+3) |n, l,m, ms) = |2,1,-1,—5). 


(b) In the presence of a spin-orbit interaction, we use the rules that j = l + 4 for l Æ 0 and j = 5 for l = 0, 
together with the fact that there are 27 + 1 values of mj, ranging in integer steps from —j to +5. 


For n = 2, the allowed values of l and s are 


l=O0orl and s=5 


For / = 0, the only allowed value of j is z, so the two possible states are: 
In, l, j, mj) = |2,0,5,-3), In, l, j mi) = |2,0, 3, +9). 
For / = 1, the allowed values of 7 are - and 3, so the possible states are: 
In, l, j my) =|2,1,5,-5), Int, j mi) = 2,1, 5,49), 
In, l, jm) =|2,1,5 73) [nl j mi) = |2,1, 3, —3) 
In, l, jmi) = |2,1,3,4+5), |n, l,j, mj} = |2,1, $, +3). 
According to relativistic quantum mechanics, energy levels depend on n and j, but not on / or mj, so each n = 2 
energy level splits into two, with four degenerate states 
nhim =|2,0,5,-5), 2,0,9,+3), 2L 2ta 
and four degenerate states 


In, ljm) = |2,1,3,-3),  |2,1,3,-5), |2,1,8,+5), 12,1, 3, +3). 


Comment In the Coulomb model of a hydrogen atom, the energy level with n = 2 comprises 8 different states 
characterized by the quantum numbers n, l, m and ms. When the spin-orbit interaction is taken into account, 
these states combine together to form new states labelled by the quantum numbers n, l, 7 and mj. It is a good 
check to note that we still have 8 different quantum states (though these are now split into two energy levels). 


Solution 7.3 The spin-orbit interaction is proportional to L - S. It is non-zero for all states with 1 Æ 0, and 
vanishes for states with / = 0. 


The Darwin term involves V?V(r), where V(r) is the Coulomb potential energy function due to the proton. It is a 
result of electrostatics that V?V (r) = 0 in regions where there is no charge (in this case, outside the proton). 
Taking the proton to be point-like, the first-order correction due to the Darwin term is only non-zero for energy 
eigenfunctions that do not vanish at r = 0. Because of the centrifugal barrier in a hydrogen atom, only the l = 0 
states meet this criterion. 


l 


Comment Each hydrogen atom energy eigenfunction involves a radial function that behaves as r’ as r — 0, so 


only the / = 0 eigenfunction is non-zero at r = 0. 
Solution 7.4 (a) We have 
T= (+8) . (+8) =E 201.626. 


SO 


(b) The perturbation due to the spin-orbit correction is of the form 


ôf = Veo(r) L - S. 
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In first-order perturbation theory, the ratio of the spin—orbit corrections in states |A) and |B} is 


EQ) _ (Al 6Eso |) 
E® (B| dH, |B) 


In evaluating a matrix element such as (A|6H5o| A), we must integrate over r, 0 and @, and also take the inner 
product with respect to electron spin. The integral over r is the same for both |A) and |B) because both these 
states have the same n and / quantum numbers, and therefore have the same dependence on r. We therefore have 


EO _ (A|L-8|A) 


E® (B|L-S|B) 
The state |A), with j = l + 4, gives 
LS\4) = (7 -£-8 ) | A) 
(Oe + §) — 10+ 1) — 5(3 + 1)) A) 


1 
2 
1 . 
=5 (+214 3-2 -1- 9) IA) 
l 
2 


so 
(A|L-$|A) = 
The state |B), with j = 1 — 4, gives 


L-as 5 (i -2°-8") |B) 


1 
= 5 (0-30 + 3) 1+ 1) - (5 + 1)) |B) 
1 
=a 4; —-V—1-3)|B) 
I+1 
so 
(B|E-8|p) = -H2 
Consequently, 
BY ıl 
A = 
EY) l 


Comments (1) On multiplying out the brackets in (L + -S): iL + 5), we have used the fact that L-S = S . È. 
This is because the orbital angular momentum operators Lz; L and L, (which act on the spatial functions) 
commute with the spin angular momentum Sas S, and S- (which act on spin kets). 


(2) In the presence of the spin—orbit interaction, the energy eigenstates of a hydrogen atom are labelled by the 
quantum numbers n, l, 7 and mj. Moreover, the spin quantum number of all states is s = 2. So each energy 


~2 a2 a2 a A 
eigenstate is an eigenfunction of J , L and S , and that is why it is useful to express L - S in terms of these 
operators. 

(3) The spin-orbit interaction increases the energy of the state with 7 = l + 5 and decreases the energy of the state 
with j =l- 3. 
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Topic 8 — Many-electron atoms 


Solution 8.1 In the central-field approximation, each electron is in a spherically-symmetric potential energy 
well. This allows the time-independent Schrödinger equation to be separated; the energy eigenfunctions can be 
written as products of a radial function R(r) and an angular function, Y (0, ¢). The angular function satisfies the 
same angular differential equation as for a hydrogen atom, and has spherical harmonic solutions Yim (0, ¢), 
labelled by orbital angular momentum and magnetic quantum numbers / and m. The radial equation is not the 
same as for a hydrogen atom, but its solutions can still be labelled by an index n (as well as by 1). 


Comment Another way of understanding this is to note that, with a spherically-symmetric potential energy 


; ‘ : shee = : ‘ : 
function, the Hamiltonian operator commutes with L and L, and this allows the energy eigenfunctions to be 


chosen so that they are also eigenfunctions of E and L z, labelled by the corresponding quantum numbers / and m. 


Solution 8.2 (a) Forn = 1, the only atomic shell is 1s (J = 0, capacity 2 electrons). 
For n = 2, the atomic shells are 2s (J = 0, capacity 2 electrons) and 2p (l = 1, capacity 6 electrons). 


For n = 3, the atomic shells are 3s (l = 0, capacity 2 electrons), 3p (l = 1, capacity 6 electrons) and 3d (l = 2 
capacity 10 electrons). 


For n = 4, the atomic shells are 4s (l = 0, capacity 2 electrons), 4p (l = 1, capacity 6 electrons) and 4d (l = 2, 
capacity 10 electrons) and 4f (l = 3, capacity 14 electrons). 


(b) The energy ordering of these shells is: 1s, 2s, 2p, 3s, 3p, 4s, 3d, ... 
Silicon has 14 electrons so the ground-state electronic configuration is 
1s? 2s? 2p° 3s? 3p”. 
Iron has 26 electrons so the ground-state electronic configuration is 
1s 2s? 2p® 3s? 3p® 4s? 3d®. 


Comment The concept of an electronic configuration, specified by the numbers of electrons with given n 
and / quantum numbers, is based on the central-field approximation. Beyond this approximation, residual 
electron—electron interactions split the ground-state configuration into various atomic terms and the spin-orbit 
interaction further splits the atomic terms into atomic levels. 


Topic 9 — Atomic terms and levels 
Solution 9.1 The possible values of L range in integer steps from |/; — l2| to lı + l2, so the possible values are 
L=2,L=s3andl=4. 
The possible values of S range in integer steps from |s; — s2| to sı + s2, so the possible values are S = 0 and 
S=], 
The corresponding atomic terms are: 
(L=2 S= 0); (b=2.5 =—1), 
(L= 3; S =0); (L=3, S = 1), 
(D=4 5S =0), (L=48 =1). 
These are labelled, respectively, by the spectroscopic symbols 1D, °D, 'F, 3F, 1G and °G. 


Comments (1) Residual electron—electron interactions produce atomic terms with definite values of L and S' for 
the whole atom. (We are implicitly assuming that the LS-coupling approximation is valid, so this scheme is not 
undermined by the spin—orbit couplings of the orbital and angular momenta of individual electrons.) 


(2) The total number of states associated with the given configuration can be calculated from the (l, s) quantum 
numbers of the two valence electrons. This gives (3 x 2) x (7 x 2) = 84 states. We can also calculate the 
number of states associated with each of the atomic terms. 1D has 5 states, 3D has 5 x 3 = 15 states, 'F has 7 
states, F has 7 x 3 = 21 states, 'G has 9 states and 3G has 9 x 3 = 27 states. It is a good check to note that 
5+ 154+ 74+ 2149-4 27 = 84. 
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Solution 9.2 The two electrons in a helium atom are identical fermions, so must be described by an 
antisymmetric total wave function. The singlet spin state is antisymmetric, so must be accompanied by a 
symmetric spatial wave function. The triplet spin states are symmetric, so must be accompanied by an 
antisymmetric spatial wave function. 


In a state described by a symmetric spatial wave function, the two electrons huddle together and experience a 
strong positive potential energy from electron—electron repulsion. In a state described by an antisymmetric 
spatial wave function, the two electrons are segregated and experience a weaker positive potential energy from 
electron—electron repulsion. 


Consequently, the singlet state, which is accompanied by a symmetric spatial wave function, has a higher energy 
than the triplet states, which are accompanied by an antisymmetric spatial wave function. 


Comments (1) The energy splitting between singlet and triplet states can be found using first-order perturbation 
theory, where the perturbation is the repulsive potential energy between the two electrons. 


(2) This is an example of the energy splitting between different atomic terms that can result from residual 
electron—electron interactions (i.e. electron—electron interactions ignored in the central-field approximation). 


Topic 10 — Molecules 


Solution 10.1 A diatomic Bez molecule has 8 electrons. Each o molecular orbital holds up to 2 electrons, so the 
ground-state configuration of a Bez molecule is 


24 ,29,25,2 
log loy 205 204- 


There are four electrons in bonding orbitals and four electrons in antibonding orbitals, so the formal bond order is 
(4 — 4)/2 =0. 


We would not expect the Bez molecule to be strongly bound. Because we are only working within the LCAO 
approximation, we can expect either a very weakly bound molecule, with a small dissociation energy, or an 
unstable molecule, but cannot be more precise than this. 


Comment As shown in Table 6.2 of Book 3, the spectroscopic binding energy of Beg is only 0.10 eV, which is 
an order-of-magnitude smaller than the spectroscopic binding energies of homonuclear diatomic molecules 
formed from atoms near Be in its row of the Periodic Table. 


Solution 10.2 (a) The Born—Oppenheimer approximation is based on the classical notion that nuclei are 
much more massive than electrons and so are more sluggish. So far as the electrons are concerned, the 
Born—Oppenheimer approximation treats the nuclei as being in fixed positions. For a given set of nuclear 
positions, the time-independent Schrödinger equation for the electrons is solved (approximately if necessary) to 
find the electronic energy eigenfunctions and eigenvalues. For each electronic energy eigenvalue, the eigenvalue is 
added to the mutual potential energy function of the nuclei to give the total static energy of the molecule. The total 
static energy is a function of the nuclear positions, and this function provides an effective potential energy function 
which appears in the time-independent Schrödinger equation associated with the motions of the nuclei. There are 
different nuclear time-independent Schrödinger equations for each electronic energy eigenvalue. 


(b) Fora given electronic energy eigenvalue in a diatomic molecule, an energy curve is a plot of the total static 
energy as a function of the internuclear distance. Assumption of the Born—Oppenheimer approximation allows us 
to calculate the electronic energy eigenvalue, and hence the total static energy, of a molecule at a series of fixed 
internuclear distances. 


Comment The Born—Oppenheimer approximation is not the same as the method of separation of variables. 
The time-independent Schrédinger equation for a molecule cannot be separated into independent equations 

for the electrons and nuclei, but the Born—Oppenheimer approximation does the next best thing. The nuclear 
coordinates appear in both the electronic and nuclear time-independent Schrédinger equations, but are given fixed 
values in the electronic time-independent Schrödinger equation. The electronic coordinates appear only in the 
electronic time-independent Schrödinger equation; in the nuclear time-independent Schrödinger equation the 
electronic degrees of freedom are replaced by an electronic energy eigenvalue (which depends on the nuclear 
positions). The coupling between the electronic and nuclear motions is revealed by the existence of different 
nuclear time-independent Schrödinger equations for different electronic energy eigenvalues. 
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Solution 10.3 (a) The total static energy is the sum of the electronic energy eigenvalue and the proton-proton 
repulsion energy, and so is given by 


O e? 1 at 1 
— 4meg \5 RFR) 
(b) The total static energy, as a function of R, provides an effective potential energy function in which the 


protons move. The equilibrium proton-proton separation is the value of R for which this effective potential energy 
function is a minimum, and is given by the equation 


dE stat (2) E e? ŚL —0 
dR 4nreg\ RE R2} ” 


which has the solution R = a. The spectroscopic dissociation energy Deg is equal to the total static energy at 
infinite separation minus the total static energy at the equilibrium proton-proton separation, and is given by 


Deg =0- Estat(a) 


o @ fle 1 
-~ 4reo \5a® a 


e2 


5rega` 


Comment The directly-observed dissociation energy is slightly smaller than the spectroscopic dissociation 
energy because the protons have zero-point vibrational energy within the effective potential energy well, provided 
by Esta( R), making the minimum energy slightly greater (i.e. less negative) than Estala). 


Solution 10.4 (a) To ensure that 7; is normalized, we must have 
1 = (pii) 
= |C? (Cos + (oil) (lis) + I4i5)) 
= |C? (ARIPI) + (oislots) + (PRIA) + (Piles) - 


Because the 1s orbitals centred on the two protons are normalized, we have 
(oilor) = (PislGis) = 1- 

Because the 1s orbitals are described by real-valued functions, we have 
(Pile) = (PRIOR): 

Hence 
1 = 2|CP (1+ (islets) 

so a suitable value of C is 
C = ——.. 

2(1 + (Gi5/¢%s)) 


(b) The function 7;(r) leads to a high probability density of finding the electron in the region between the 
protons. This partially screens one proton from the other and promotes bonding. This mechanism does not apply 
for the function 72(r), which has a node midway between the protons. 


Comment Note that atomic orbitals centred on different protons are not orthogonal, so it would be wrong to take 


C=1/,/2. 
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Topic 11 — Solids 


Solution 11.1 Bloch’s theorem states that, in an infinite crystal, the energy eigenfunctions of electrons can be 
written in the form 


v(t) = e™Tux(r), 
where k is a real, constant vector and u,(r) is a periodic function, which may depend on k and has the periodicity 
of the lattice: 


ux(r +R) = ux(r) for any lattice vector R. 


Lattice points are separated by a lattice vector R, and Bloch’s theorem gives 
lV +R)? = juf +R)? = [u(e) |? = |W). 


Comment Each Bloch energy eigenfunction is the product of a periodic function (which is the same at each 
lattice point) and a plane wave function (which varies from one lattice point to another, unless k = 0). A typical 
electronic energy eigenfunction 7 does not have the periodicity of the lattice, but the probability density |7)|? does 
have this periodicity because the plane wave factor has unit modulus. 


Solution 11.2 (a) For any lattice vector R,, replacing r by r + R; throughout u(r) gives 
i 


= 37 o(r— (R; — Rye RBs, 


Here, R; — R; is another lattice vector, Rj. Because the sum is over all lattice vectors in an infinite crystal, we 
obtain 


u(r + Ry) = X g(r — Rye = u(r). 
l 


Since u(r) = u(r + R4) for any lattice vector R;, we can say that u(r) has the periodicity of the lattice. 


(b) Multiplying both sides of the expression for u(r) by e€ 


u(r)e — eik-r » e(r = Rj) 


gives 


= > olr = R, je Ri 


= V(r), 
so we conclude that y(r) is the product of a plane wave e 
and therefore satisfies Bloch’s theorem. 


ikr and a function u(r) with the periodicity of the lattice, 


Comment The function u(r)e'* is used in the tight-binding approximation for the electronic energy levels in a 
crystalline solid. It is used as the trial function in the variational method, even though there are no adjustable 
parameters in this case. 


Solution 11.3 For a hydrogen atom, the Rydberg energy is given by 
2 2 
e H 
ER = | — | > = 13.6eV. 
x — ) 2h? i 

The corresponding result for a donor atom in germanium is obtained by replacing £o by 15.8£0, and jz by the 
effective electron mass mz = 0.12 me, so the energy that must be supplied is 

0.12 


D” 13.6eV = 6.5 x 107° eV. 


Comment The typical thermal energy at room temperature is 0.025 eV, which is much larger than the ionization 
energy we have calculated. So, at room temperature, practically all pairs of donor atoms and electrons are ionized. 
At much lower temperatures the donor atoms and electrons form bound states, leading to a drop in the number of 
electrons in the conduction band and a drop in the electrical conductivity. 
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Solution 11.4 A metal has a partly-filled energy band and a fixed number of current carriers. Its electrical 
conductivity is limited by deviations from a perfectly periodic crystal lattice. At normal temperatures a significant 
contribution from these deviations arises from vibrations of the crystal lattice, which increase with increasing 
temperature. The electrical conductivity of a metal therefore decreases with increasing temperature. 


In a pure semiconductor at absolute zero, the valence band is full and the conduction band is empty, but the gap 
between the valence and conduction bands is small enough that, at higher temperatures, a significant number of 
electrons can be excited into the conduction band, leaving holes in the valence band. Conduction then takes place 
due to both electrons in the conduction band and holes in the valence band. The number of these electrons and 
holes increases with increasing temperature. This generally outweighs any effect due to an increase in lattice 
vibrations, so the electrical conductivity of a typical pure semiconductor increases with increasing temperature. 


Comment In a doped semiconductor, there may be a range of temperatures over which almost all the impurities 
have donated electrons to the conduction band or produced holes in the valence band, but for which hardly any 
electrons have been excited from the valence band to the conduction band of the host semiconductor. In such a 
case, the electrical conductivity of the semiconductor may decrease with increasing temperature, as for a metal. 


Topic 12 — Interaction of light and matter 


Solution 12.1 (a) Absorption and stimulated emission have transition probabilities that can be predicted 
directly from first-order time-dependent perturbation theory. 


(b) Quantum field theory is required for a fundamental understanding of spontaneous emission. In quantum field 
theory, the vacuum is not quiescent, but is a seething soup of virtual particles which can disturb an atom and shake 
it out of an excited stationary state. 


Comment First-order time-dependent perturbation theory allows us to calculate a transition probability 
due to the presence of an external time-dependent perturbation. In absorption and stimulated emission the 
time-dependent perturbation is due to an externally-applied electromagnetic wave. 


Solution 12.2 According to the results of time-dependent perturbation theory, the transition probability from an 
initial state 7); to a final state wf is proportional to 


[dbp] VCE) lal? œ Isl cos 8 lyi)? 
for light polarized along the z-direction. 


For a hydrogen atom, the initial and final states are of the form 
vilr, 0,0) = Rng; (7) O1,m, (Ae? 
prit 4) = Ray (7) O1pm, (Ae? 


The transition probability therefore depends on the matrix element 


27 fT oo 
(jl cosas) = f [| W7(r, 0, p) 7: cos 0 yilr, 0, p) r° sin dr d0 do 


oO 


2T T 
a eilm—mp)8 ag | GF m, (0)Oim; (0) cos sin 0 d0 f * (7) Raat (r) dr. 
0 0 0 


For m; # my, the integral over ¢ gives 


; p=2r 
2m i(mi=mys)o 
f eilmi—™5)$ qdo = — = 0, 


0 My TF) 2g 


so the radiative transition will not take place (in first-order time-dependent perturbation theory) unless ms = mi. 
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Comment For radiative transitions induced by electromagnetic radiation polarized in the z-direction, there is 
a selection rule forbidding transitions with mf # m. Note that this selection rule applies only to radiative 
transitions (it does not apply to transitions induced by atomic collisions) and it is restricted to transitions induced 
by electromagnetic radiation polarized in the z-direction. Moreover, it is only valid in first-order perturbation 
theory. Transitions with ms Æ m; may still occur in this case, but they are strongly inhibited and may not be 
noticeable compared to transitions with m+ = m. If they occur at all, they will produce very faint spectral lines. 


Solution 12.3 The solution follows the same lines as Exercise 12.2, but with 
V(t) = f(t) rsin 8 cos ¢. 


In this case, the relevant matrix element is 


20 T lo) 
(jl sin cos 6) = | [| We(r, 0, d) r sin 8 cos 6 Y(r, 0, $) r° sin 6 dr dO do 


27 T love) 
- I ellm-m7)9 cos ddd / GF sn; (0)Oim; (0) sin? 6 d0 f * y (r)Run (r) dr. 


For m; = mp, the integral over ¢ gives 


ma cos ọ dọ = [sin J au = 0, 
0 ? 


so the radiative transition will not take place in this case (in first-order time-dependent perturbation theory). 


Comment For electromagnetic radiation polarized in the x-direction, the selection rule forbids transitions 

with m; = my (cf. Exercise 12.2). Again, this selection rule applies only for radiative transitions and for a 
particular polarization. It is valid within first-order perturbation theory, and so may not completely forbid very rare 
transitions. 


Solution 12.4 (a) The three types of radiative transition are as follows: 


(i) In absorption, a system absorbs a photon and makes a transition from a quantum state of lower energy to a 
quantum state of higher energy. 


(ii) In stimulated emission, a system is stimulated by an incident photon to make a transition from a quantum 
state of higher energy to a quantum state of lower energy. In the process, a photon is emitted with the same 
characteristics as the incident photon. 


(iii) In spontaneous emission, a system emits a photon and makes a transition from a quantum state of higher 
energy to a quantum state of lower energy, irrespective of any incident radiation. 


(b) The term Bı2u(w21)Nı corresponds to the rate of transfer of atoms from quantum state 1 to quantum state 2 
due to absorption. 


The term B2)u(w21) N2 corresponds to the rate of transfer of atoms from quantum state 2 to quantum state 1 due 
to stimulated emission. 


The term A21 Nə corresponds to the rate of transfer of atoms from quantum state 2 to quantum state 1 due to 
spontaneous emission. 


Comments (1) Note that the first index in B12, B21 and Ag refers to the initial state and the second index refers 
to the final state. 


(2) Using the Boltzmann distribution law and the Planck distribution law, we can show that the two B-coefficients 
are equal to one and obtain an expression for the A-coefficient in terms of the B-coefficients. The B-coefficients 
are found from first-order time-dependent perturbation theory, and the A-coefficient found from them. This gives 
us an indirect way of predicting the rate of spontaneous transitions from an upper quantum state to a lower 
quantum state. 
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